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Abstract. We present a model for the aperiodic variablity of accreting 
Black Hole Candidates (BHC) along with model light curves. According 
to the model this variability is the combined outcome of random (Poisson) 
injection of soft photons near the center of an extended inhomogeneous 
distribution of hot electrons (similar to those advocated by the ADAF 
or ADIOS flows) and the stochastic nature of Compton scattering which 
converts these soft photons into the observed high energy radiation. Thus, 
the timing properties (PSD, lags, coherence) of the BHC light curves reflect, 
to a large extent, the properties of the scattering medium (which in this 
approximation acts as a combination of a linear amplifier/filter) and they 
can be used to probe its structure, most notably the density profile of the 
scattering medium. The model accounts well for the observed PSDs and 
lags and also the reduction in the RMS variability and the increase in 
the characteristic PSD frequencies with increasing source luminosity. The 
electron density profiles obtained to date are consistent mainly with those 
of ADIOS but also with pure ADAF flows. 



1. Introduction 

The study of the physics of accretion powered sources, whether on galac- 
tic (X-ray binaries) or extragalactic systems (AGN), involves length scales 
much too small to be resolved by current technology. As such, this study 
is conducted mainly through the theoretical interpretation of their spec- 
tral and temporal properties. Until recently studies of this class of objects 
focused, for technical mainly reasons, on their X-ray spectra. These are 
generally fit very well with those of Comptonization of soft photons by 
hot (Te < 10^ K) electrons, a process explored in great depth over the past 
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twenty or so years ([11]; R. Sunyaev this volume). Because the electron tem- 
peratures of matter accreting onto a black hole are expected to be similar 
to those necessary to produce the observed spectra, it has been considered 
that detailed spectral fits of these sources would lead to insights on the 
dynamics of accretion onto the compact object. 

However, the determination of accretion dynamics requires the knowl- 
edge of the density and size of the emitting region, neither of which is pro- 
vided by radiative transfer and spectral fitting considerations (the equations 
of radiative transfer involve the optical depth as the independent variable) . 
Indeed, as shown explicitly in [6] and [5], plasmas of very different radial 
extent and density profiles can yield identical Comptonization spectra. The 
degeneracy of this situation can lifted with the additional information pro- 
vided by timing observations. 

The timing properties of BHC, however, suggest length scales inconsis- 
tent with the prevailing notion that the observed X-rays are emitted from 
a region of size a few Schwarzschild radii, Rs' The power spectra (PSD) 
of BHC exhibit most of their power at scales ~ 1 s, far removed from the 
characteristic time scales associated with the dynamics in the vicinity of 
the black hole horizon, Rs/c ~ 10"'^ s. Until recently rather little attention 
has been paid to this time scale discrepancy, generally attributed to the 
(unknown) mechanism "fueling" the black hole, presumably operating at 
much larger radii; rather, more attention was paid to the fiicker noise-type 
(oc f~^) PSDs of this class of sources. Novel insight into the variability 
properties were introduced by [8], who showed that both the magnitude 
and the Fourier period dependence of the lags in the X-ray light curves at 
two different energies was inconsistent with Comptonization by a plasma of 
size ~ a fewRs; the magnitude of time lags (~ 0.1 s) suggested an emission 
region of much larger size, precluding thus an explanation of the observed 
PSDs due to a modulation of the accretion rate onto the black hole. 

These discrepancies between the expected and the observed variability 
of GBHC led [6], [3], [5] and [7] to propose that, contrary to the prevailing 
notions, the size of the scattering region responsible for the X-ray emission 
is not i? ~ 3 — 10 i^s but rather R ^ lO^Rs, as implied by the PSDs and 
lag observations. Furthermore, the scattering medium (corona) is inhomo- 
geneous, with the electron density following the law, 



where r is the radial distance from the center of the corona (assumed to be 
spherical) and ri, r2 are its inner and outer radii respectively. The index 
p > is a free parameter whose value depends on the specific dynamical 
model that determines the electron density. For example, the ADAF of [10] 
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Figure 1. The response function g{t) of a corona with ri — 10^^ light s, r2 — 10 L s, 
p = 1, tt = 1- The sohd hne is a fit to the MC data (dashed hne). The dotted hne is the 
Fourier spectrum |G(a;)p of g{t). 



suggest p = 3/2 and Tg ^ 10^ at radii as large as r ~ {mp/me)Rs, while 
models combining inflow and outflow [2], [1] (ADIOS) allow in addition 
values p ~ 1. As pointed out in [6], [5] most of the data analyzed to date 
suggest p = 1 with p = 3/2 also acceptable in certain cases. 

It is intuitively obvious that scattering in the extended configuration 
of Eq. (1) produces time lags over a range of Fourier periods similar to 
the range of radii span by the hot corona: Scattering at a given radius 
R increases the X-ray energy and introduces a lag At ^ i?/c between the 
scattered and unscattered photons, the lag appearing at a Fourier period 
P = R/c ([5]; N. Kylafis these proceedings). To compare this to the lags 
given in [8] (which are the average over all such photon pairs scattered in a 
given decade in R, as a function of R) one has to multiply At by the proba- 
bility of scattering at a given radius, ViR) — t{R), with t{R) the scattering 
depth over the radius R. For the density profile of Eq. (1), t{R) oc R~p~^^ 
and since the Fourier period P oc R, (At) oc R~P'^'^ oc P~^'+^. Monte Carlo 
simulations and analytical considerations have shown that these arguments 
are essentially correct [5]. In addition, these simulations showed that the 
configuration of Eq. (1) produces light curves of high coherence over the 
range [r2/c, ri/c] of the Fourier period, provided that the soft photons are 
injected near its center. 



2. The Model Light Curves 

Based on these considerations one can easily produce model light curves, 
whose properties in the time or the Fourier domain can then be compared 
to observations for consistency. The high coherence of the observations [12] 
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alpha=0.5, f=3 — p 
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Figure 2. The light curve F{t) of Eq.(2) for the g{t) of Figure 2 and / = 3. 

suggest injection of the soft photons at r << r2. The response function, g{t), 
of a corona with density profile given by Eq. (1) to an instantaneous release 
of soft photons at r < ri has been computed by a Monte Carlo simulation 
[6], [7]. The result of a particular case is shown in Figure 1 along with an 
analytic fit. As pointed out in [5] these functions can be approximated well 
by a Gamma function distribution i.e. g{t) oc t"^^ e^*/*^ (0 < a < 1, /3 > 0). 

Assuming linearity, i.e. that Tg is not affected by the photon flux, model 
light curves can be produced by an incoherent, random (Poisson) injection 
of shots of the form g{t). The prescription for such a light curve is 

F{t)=J2Qi0(t-J2ii]9{t-ti) . (2) 

i=l \ i=l J 

Qi are the shot amplitudes, assumed constand, 9{t — ti) is the Heaviside 
function and ti is a collection of Poisson distributed time intervals obtained 
from the expression ti = —f ■ to ■ logRi] Ri is a random number uniformly 
distributed between and 1, and / is a real number, indicating the mean 
time between shots in terms of their rise time to. Figure 2 shows the light 
curve obtained using the above prescription with g{t) as given in Figure 1 
and / = 3, while Figure 3 the light curve from an identical sequence of Ri's 
and / = 10. 

These light curves look similar to those of BHC sources in their low, 
hard state [9]. Due to the Poisson nature of tj's and the identical shape 
of the shots that make up these light curves, their PSDs are those of an 
individual shot, also shown in Figure 1 (dotted line). The form of the PSD is 
also quite similar to those of BHC (see [9]), with the high and low frequency 
breaks associated respectively with the inner ri and outer r2 edges of the 
scattering corona. 
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Figure 3. Same as Figure 2, with / = 10. 



One can easily produce, in addition, light curves for different photon 
energies, considering that for photons of a higher energy, g{t) is approxi- 
mately of the same form with only a small increase in the parameters a 
and (5 [5]. While the resulting light curves are visually identical to those 
of lower energies, their difference is nonetheless easily seen in the phases of 
the corresponding FFT. We have generated this way model light curves ap- 
propriate to three different energies and then computed the corresponding 
pairs of phase lags (a)-(b), (a)-(c) as a function of the Fourier frequency 
measured in Hz; the results are shown in Figure 4. The phase lags as pre- 
sented in this figure are very similar to those of the sources analyzed in [9] 
and in particular the X-ray transient GRO J0422+32 [4]. 

A correlation between the RMS variability and the value of / is the most 
prominent feature of Figures 2 and 3. Because, in the linear regime, smaller 
/ is equivalent to higher luminosity (more photons per unit time), this cor- 
relation provides a natural account of the observed anticorrelation between 
the BHCs luminosity (in their hard spectral states) and their RMS variabil- 
ity. Furthermore, consideration of this model within its natural dynamic 
framework, namely that of ADAF, provides a direct correspondence be- 
tween the dynamical and the timing properties of these systems: In ADAF 
all scales (and also ri and r2) scale proportionally to VffTcooi with Vff the 
free-fall velocity and TcooI the local cooling time. For r^^oo/ inversely propor- 
tional to the local denisty (as is the case with ADAF), all length scales 
(and therefore frequencies) associated with the corona should depend on 
rn/niEdd-, decreasing with increasing value of rh/rriEM'i hence the corre- 
sponding frequencies should increase with increasing luminosity, a behav- 
ior which has apparently been observed in most accreting sources (whether 
neutron stars or black holes). 
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Figure 4- The phase lags of two sets of model light curves with the following parameters 
(a) a = 0.5, P = 16 sec, (b) a = 0.55, /3 = 16 sec, (c) a = 0.6, /? = 16 sec. The two curves 
correspond to the lags between: (a) - (b) (solid line) and (a) - (c) (dotted line). 



3. Conclusions 

(a) The model of the extended inhomogeneoiis corona provides in a simple, 
well understood fashion model light curves for BHC which have morphology, 
PSDs and phase (or time) lags very similar to those observed, (b) The 
dependence of the phase lags on Fourier period allows the determination 
of the density profile of the corona; the data are consistent both with p = 
3/2 (ADAF) and p = 1 (ADIOS) flows, (c) The dependence of the RMS 
variability and the PSD break frequencies of BHC on their luminosity are 
also in general agreement with these dynamical considerations and provide 
a well defined framework within which these ideas could be tested through 
more detailed comparison with observations. 
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